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The paper by Sato and Kobayashi in 1977 studied the cosmological effects of a massive 
neutrino and obtained constraints on its properties. This paper initiated many studies to 
use cosmology as a laboratory of particle physics or to use particle physics to explore the 
very early universe. 

§1. Introduction 

In 1977 Kobayashi and one of the present author (K.S.) studied cosmological 
effects of a massive neutrino and published the paper entitled "Cosmological Con- 
straints on the Mass and the Number of Heavy Lepton Neutrinos'' Hereafter we 
call it Sato-Kobayashi 1977. This study was motivated by the discovery of anoma- 
lous lepton production in e + -e~ annihilation (e + + e~ — > e^ + ^ + missing energy) 
in 1975, 2 ) which strongly suggested the existence of a heavy lepton having a mass 
in the range 1.6 to 2.0 GeV. [ Later, the existence of new lepton was confirmed and 
named tau (r). ] 

Assuming the heavy lepton in the third family, the Sato-Kobayashi 1977 ad- 
dressed two issues; one was to constrain the mass of the associated neutrino and the 
other was limiting the number of massive neutrinos. The latter issue had been partly 
investigated by Steigman et al. 3 ' for massless neutrinos and Kobayashi-Sato 1977 
gave a more generic limit which applies to massive neutrinos. As for the first issue 
the authors were first to study the cosmological effects of a massive neutrino and ob- 
tained constraints on its mass and lifetime. In deriving the cosmological constraints, 
the effects on the cosmic matter density, cosmic background radiation and big bang 
nucleosynthesis were considered. Today this cosmological approach is the standard 
way to reveal properties of particles predicted in new models of particle physics or 
to explore the very early universe by using new particle theories. The research field 
is now called particle cosmology. Therefore, Sato-Kobayashi 1977 is regarded as a 
historic paper that took the initiative in development of particle cosmology. 

In this paper we review Sato-Kobayashi 1977 and the subsequent progress. The 
paper is organized as follows. In section [2] we review Sato-Kobayashi 1977. In 
section [31 the subsequent progress is shown, and we conclude in section 01 
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§2. Review of Sato-Kobayashi 1977 

2.1. Lifetime of the massive neutrino 

In the 1977 paper, 1 ) a neutrino, which is associated with the heavy lepton ( now 
called 'tau') suggested in the e + -e~ collider experiment, was considered. If it has a 
finite mass and qurk-like mixings with other neutrinos (f e ,z/„), the heavy neutrino 
is unstable and can decay into a lighter neutrino via 

»h -> ^ + 7- (2-1) 

The neutrino with mass m Uh > 2m e can also decay into an electron-positron pair as 

v h ^u e + e + + e~. (2-2) 

With use of the mass eigenstates v% (i = 1,2,3) the weak interaction eigenstate 
of the heavy neutrino is written as 



J2 U hM> ( 2 ' 3 ) 



where is the neutrino mixing matrix. In Ref.l), the lifetimes for the decay 
processes (|2-ip and (|2-2p were estimated in the framework of Weinberg-Salam theory. 
For decay mode (|2-1[) 

128tt 3 16 TT TTif ml \ ~ 2 

* 2xl *Hw)~ 6 v^~ 2 > (24) 

where m Uh is the mass of the heavy neutrino, rrii is the charged lepton mass, a is the 
fine structure constant, Gf is the Fermi constant and M\y is the W boson mass. In 
the second line of Eq. (|2-4I) we have used the measured r lepton mass, but in 1977, 
the heavy lepton (= r) was not confirmed, so the upper bound on the lifetime was 
obtained from unitarity of the mixing matrix, which leads to the following inequality: 



i'<i j j 

1 



2\2 



<-max[(m?-m?,) 2 ]. (2-5) 
4 jd , 



Using this relation and rrij/Mw < 1 5 the lower bound on the lifetime was obtained 

as 

-2 



128vr 3 64 /max,- [774] \ " / m \-5 
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On the other hand, the lifetime for decay into e + e is given by 

( J- +^ ^ 1927rS ITT TT* 1—2 

r(i/ h ^^ e + e^ + e ) = — —\Ui h U le \ 

- 3Xl04seC (Sv)"l^"" 2 ' ^ 

for > 2m e . The mixing is restricted with use of fi-e universality. The pion 
decay process tt — ► e+u^, if exists, is enhanced because the amplitude is proportional 
to m Uh . So in Ref. 1) requiring that the contribution of tt — > e + v^ to the total decay 
width of tt — > e + f should be less than 5%, the limit on the U\h was estimated as 

\U lh \ 2 (m Uh /m e ) 2 < 0.05, (2-8) 

which led to the lower limit to the lifetime, 

r{y h - v e + e + + e") > 2 X 10 8 sec (^) . (2-9) 

Eqs. (l2-6p and ()2-9[) show that the massive neutrino decays with cosmological 
time scale and the decay into an electron-positron pair, if kinematically possible, 
may be faster than that into a photon. 

2.2. Cosmological constraints 

Since the massive neutrino is long-lived (r > 1 sec) , it can affect the thermal his- 
tory of the universe. In Sato-Kobayashi 1977 1 ) the constraints on mass and lifetime 
were obtained by considering the following cosmological effects: 



1. When the massive neutrino is stable or has a longer lifetime than the present 
age of the universe, its present mass density may exceed the observed upper 
limit. 

2. Photons and electrons produced in the decay may cause a spectral distortion 
of the cosmic microwave background (CMB) radiation, or if the decay takes 
place after recombination, the emitted photons may be directly observe as the 
background radiation. 

3. The massive neutrino can affect the big bang nucleosynthesis (BBN) by speed- 
up of the cosmic expansion and/or entropy production. 



First let us consider the cosmic density of the heavy neutrino. The neutrinos are 
in thermal equilibrium via weak interaction {y^ + i>h e~ + e~) in the very early 
universe. The annihilation cross section is 

f GlT 2 m Vh < T 

m~ * h : T , (2-io) 

{ G 2 F mi h m Uh > T 

where v is the relative velocity and (• • • ) means thermal average. When the interac- 
tion rate r = (o~v)n I/ becomes less than the cosmic expansion rate H ~ T 2 /M p (M p : 
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Planck mass ~ 2.4 x 10 18 GeV), the neutrino decouples from the thermal bath. Af- 
ter decoupling, the neutrino does not interact with other particles and its comoving 
density freezes out. The freezing temperature is determined from r ~ H, leading to 
Tf ~ max[l, (m^/lOMeV)] MeV. Then the number ratio of the massive neutrino to 
the photon is estimated as 

\ Q /III T\/T^r\ 3 I 

(2-11) 

Notice that the neutrino number density decreases for m Vh > 10 MeV because such 
heavy neutrino becomes non-relativistic at decoupling and its number density suffers 
from the Boltzmann suppression (~ exp(— m/T)). 

2.2.1. Stable Neurino 

If the massive neutrino is stable or has a longer lifetime than the age of the 
universe, the neutrino density should be less than the upper limit of the present 
density of the universe. The present mass density of the heavy neutrino is given by 

3 f /10MeV\ 3>l 



Pu h = —g s m Uh n li0 mmy., — — J j , (2-12) 

where g s is the statistical weight of the massive neutrinoEH In term of the density 
parameter i? = p/p c , where p c is the critical density (= 1.054 x 10 4 /i 2 eV/cm 3 ) and 
h is the Hubble constant in units of lOOkm/s/Mpc, the above equation is rewritten 
as 

, 2 / rn Uh \ f / 10 MeV\ 3 l 



In Ref.l), from observations of the deceleration parameter [qo = (da /a 2 )] and the 
age of the universe at that time, f2h 2 < 1.45 was adopted as the upper limit of the 
density parameter, which resulted in 

m Uh < 70 eV or m Uh > 3.7 GeV. (2-14) 

As mentioned before, since the Boltzmann suppression reduces the number density 
of the heavy neutrino, the neutrino with large mass vn Vh > 3 GeV is allowed as well 
as the light neutrino with mass less than O(10) eV. Unfortunately, the lower bound 
on the heavy neutrino mass is often called "Le e- W einberg limit" because Lee and 



Weinberg 5 ^ obtained the limit at the same timer*) 



*' The statistical weight g s is always 1 for a Majorana neutrino. As for a Dirac neutrino g s is 
1 — 2 depending on its mass. Because of chirality of the weak interaction, the rate of the interaction 
of a right-handed neutrino is suppressed by a factor (m v /E) 2 where E is the relevant energy. Thus, 
neutrinos with mass less than 2 keV are never in thermal equilibrium in the early universe and 
g s = 1 for that case. 4 ' However, in Ref.l) g a = 2 was assumed for any massive neutrinos. 

**' The both papers 1 '' 5 ' were submitted in May 1977. We should also note that the lower bound 
on the heavy neutrino mass was also obtained by P.Hut, 6 ' Dicus, Kolb and Teplitz 7 ' and Vysotskii, 
Zel'dovich and Dolgov. 8 ' 
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2.2.2. Unstable Neutrino 

Next, we consider the constraint on the unstable massive neutrino which de- 
cays into a photon or an electron-positron pair. If the neutrino decays sufficiently 
early t < 10 10 sec, the emitted photons are quickly thermalized through Compton 
scattering and bremsstrahlung. However, when decay takes place at later epochs, 
the emitted photons may partially thermalized via Compton scattering but cannot 
form the Planck distribution because photon number changing processes such as 
bremsstrahlung are inefficient. Thus, the CMB spectral shape is distorted, which 
leads to a constraint on the lifetime r Vh of the massive neutrino. 

When photons are emitted in the neutrino decay, the ratio of the extra photon 
density Z\p 7 to the CMB density p 7 at the decay time is estimated as 

Ap-rM _ {3/ll)g 8 {m v j2)n^F m Uh 



P t (t„J 2.7T d n n 



0.05 5s -^F, (2-15) 



where is the cosmic temperature at t = T Uh and F = min[l, (10 MeV/m^) 3 ]. 
With relation between the cosmic temperature and time, t = 10 12 sec(T/eV) -3 / 2 , 
and Ap 1 /p 1 < 0.1, the constraint was obtained 1 -* as 

' 10 12 sec(^)- 3/2 m Uh <20eV 

10 10 sec 20 eV < m„ h < 2 GeV . (2-16) 



10 9 sec(^) 3 m„ h >2GeV 

The above constraint obtained for decay into photons also applies to the case where 
electron-positron pairs are produced in the decay, i.e, for m Uh > 2m e because the 
energetic electrons (positrons) scatter off the CMB photons to produce energetic 
photons or the positrons annihilate into photons. 

When the decay takes place after recombination the emitted photons freely 
stream and may be observed as the background radiation after suffering redshift. 
This extra background radiation should not exceed the observed one, from which 
Sato-Kobayashi 1977 obtained the constraint on the lifetime of the massive neutrino 
as 

T Vh > 10 18 sec (2-17) 

Finally, let us show how Sato-Kobayashi 1977 derived the constraint on the 
lifetime and mass of the massive neutrino from consideration of nucleosynthesis in 
the early universe. Two effects both of which increase the 4 He abundance were 
studied in Ref.l). 

The first effect is entropy production by the neutrino decay. It is well known that 
the abundances of light elements synthesized in the standard BBN only depend on 
the baryon density of the universe which is parametrized by "photon-baryon ratio" , 

VB = — = 2.7 x 10- 8 n B h 2 . (2-18) 

n 7 

In Ref.l) the different parameter = ps/Tg (Tg = T/10 9 K) was used but here we 
have adopted the standard parameter used in recent literature, hjy is related to r\B 
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as 

h N = 3.4 x 10- 6 (77 B /10~ 10 )g/cm 3 . (2-19) 

When the neutrino decay produces many photons, the baryon-photon ratio at BBN 
epoch (?7bbn) is larger than the present one and is roughly given by 

— + (2 ' 20) 

where Ap^/p^ is given by (|2- 15|) . Using the lower bound on the present baryon den- 
sity (i?B > 0.05), which predicts the lowest 4 He abundance because the theoretical 
prediction of 4 He abundance is an increasing function of r?bbn> Pbbn is estimated as 

.-10 A , A Pl( T Vh, 



%bn < 3.5 x 10- iU 1 + r y " h/ . (2-21) 

The other effect is speed-up of the universe due to the massive neutrino. As 
compared with massless case, the massive neutrino can have a larger energy density, 
which increases the expansion rate of the universe and leads to earlier freeze-out of 
neutron-proton exchange interactions (e.g. n + v e <-> p + e~). Since the neutron- 
proton ratio at the freeze-out temperature Tf is given by n/p ~ exp(— 1.29MeV/Tj), 
the earlier freeze-out (higher Tf) results in more abundance of n and hence more 
production of 4 He. The change of the cosmic expansion rate is represented by £ 
defined as 

^(^f- < 2 ' 22 > 

\Pstd{T)J 

where p st d is the cosmic density at temperature T for the standard case without 
massive neutrinos. The densities p s td and p are given by 



PmA - |j fy/'' • ^ <' j ■ P - /'s,a + (in,, + 3.157;,)/,,,,. (2-23! 



where T u and iVgtd are the neutrino temperature and the number of massless neu- 
trino species. As mentioned before, in 1977 the existence of the tau lepton was not 
confirmed, thus iV s td = 2 was taken as the standard value. Since the speed-up ef- 
fect is important when the neutron-proton ratio freezes out (Tf ~ 1 MeV.), we can 
assume T > m e and hence T v = T. Thus £ is written as 

f = l + 0.062 5s &L + 3.15) . (2-24) 

On the other hand the theoretical prediction for the 4 He abundance was calculated 
numerically by Wagoner 9 ) and he gave the empirical formula, 

Y = 0.264 + 0.0195 log(?? bbn /10~ 10 ) + 0.380 log f , (2-25) 

where Y is the mass fraction of 4 He, Y = P^uc/pb- Using Eqs. (|2-2ip . (|2-24j) and 
(|2-25p together with observational upperbound Y < 0.29, the following constraint 
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on the mass and lifetime was obtained: 1 ) 

A 1 ^(^r 2 K<lMeV) 
T u h << 4 (2-26) 

{ 1 sec (^) {m Vh > 1 MeV) 

The cosmological constraints obtained in Ref.l) were so stringent that a large 
region in the parameter space (T nUh ,m Uh ) was excluded as shown in Fig.4 of Ref.l). 
Furthermore, together with the lower limits on the lifetime of the massive neutrino 
estimated in the Weinberg-Salam theory (|2-6p and (|2-9p . the neutrino in the mass 
range, 

70 eV < m Vh < 23 MeV (2-27) 

was found to be excluded. 

2.3. Number of Massive Neutrinos 

Furthermore, Sato-Kobayashi 1977 discussed the limit on the number of massive 
neutrinos. As seen in the previous subsection, the masses of neutrinos should be less 
than 70 eV or more than 23 MeV. So we can assume that N± types of neutrinos have 
masses less than 70 eV and N2 types of neutrinos have masses larger than 23 MeV. 

The constraint on N% is obtained by considering the speed-up effect in BBN. In 
this case, assuming massless v e and 1^, Eq. (l2-24j) is rewritten as 

f = 1 + 0.196 g s N v (2-28) 

Using Eq. (|2-25p and Y < 0.29 the constraint is obtained as 

N^SgJ 1 . (2-29) 

As for N2, suppose that v m has the longest lifetime among neutrinos with mass 
larger than 23 MeV. Then the lifetime is bounded as 

(2-30) 

where vm 2 ls the neutrino with the largest mass. Using /i-e universality Sato- 
Kobayashi 1977 obtained 

r Vm > 2.2 x 10 6 N 2 (m UN2 /MeV)- 3 sec. (2-31) 

In Ref.l), applying the BBN constraint, N2 was limited as 

m VN 9 X 7 




N ' 2 ~ I 23lfeV ' ' < 2 ' 32 ' 



§3. Progress after 1977 



Since Sato and Kobayashi first obtained the constraint on the massive neutrino 
in 1977, the constraints have been improved according to development of cosmology, 
laboratory experiments and observations. 
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3.1. Stable Neutrino 

The neutrino, if stable, cannot have a mass larger than about 70 eV from the 
cosmic density constraint. In other words, the neutrino can be dark matter of the 
universe if its mass is O(10) eV. In fact, the neutrino with 0(10) eV attracted much 
interest in cosmology in 1980's. In addition, at that time it was reported that a 
Russian group 10 ) discovered the electric neutrino mass in the tritium beta decay 
( 3 H ^ 3 He + e~ + v e ) experiment [ which was later denied by other experiments]. 
However, it was soon found that a large velocity dispersion of the neutrino erases the 
small scale density fluctuations by free-streaming. So if the neutrino is dark matter, 
very large structure such as clusters of galaxies are formed first and smaller structure 
like galaxies are formed later. This "top-down" scenario was studied extensively 
by iV-body simulations and it was concluded that the neutrino cannot be a main 
component of the dark matter. u ) 

Here, before discussing the recent cosmological bound on the neutrino mass, we 
briefly describe the experimental progress towards measurement of neutrino masses. 
The experimental evidence for massive neutrinos was first obtained by the Super- 
Kamiokande(SK) experiment 12 ) which measured atmospheric neutrinos and discov- 
ered the neutrino oscillation between fj,- and r-neutrinos. The experiment suggests 

\ml T - ml J ~ 3 x 10~ 3 eV 2 . (3-1) 

Moreover, the solar neutrino experiments 13 )' 14 ) and the reactor experiment 15 ) also 
see the oscillation between e- and fx- neutrinos and obtained 

|m 2 M -m 2 J~7xl0- 5 eV 2 . (3-2) 

Notice that the oscillation experiments can measure only the differences of neutrino 
mass squares and do not give the absolute values. The most severe mass limit is still 
obtained by the tritium beta decay experiment. The present best limit is 16 ) 

m„ e < 2 eV. (3-3) 

As for the other neutrinos, the present mass bounds are 16 ) 

m„„ < 0.19 MeV (3-4) 
m„ T < 18.2 MeV (3-5) 

Let us return to cosmological consideration. Even if the massive neutrino is a 
minor component of the dark matter, its existence can give a significant effect on the 
density fluctuations on small scales. Fukugita, Liu and Sugiyama, 17 ) comparing the 
predicted density fluctuations with that inferred from the cluster abundance, derived 
the neutrino mass limit as 

^77^<3eV, (3-6) 

where the sum is taken over the three species of neutrinos. The result of the oscil- 
lation experiments strongly suggest that the three neutrino masses are degenerate if 
their masses are larger than 0.1 eV, i.e. m Ue ~ m v ~ m„ r . 
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The neutrino mass limit was significantly improved by Wilkinson Microwave 
Anisotropy Probe (WMAP) observation which obtained the full sky map of the CMB 
temperature and thereby determined the cosmological parameters precisely. The 
neutrinos with mass > O(0. 1) eV become non-relativistic before or during recom- 
bination and they change the epoch of matter-radiation equality, which alters the 
amplitude of the temperature fluctuations through change of the gravitational po- 
tential from radiation dominated era to matter dominated one. The recent WMAP 
five year data alone gives the upper bound on the neutrino mass as 19 ) 

J2 m "i < L7 eV (95%CL), (3-7) 

and a more stringent upper bound is obtained if the other observational data (baryon 
acoustic oscillation and supernovae ) are added , 19 ) 

J2m Ui < 0.61 eV (95%CL). (3-8) 

3.2. Unstable Neutrino 

In the 1977 paper, the cosmological constraints on the radiatively decaying neu- 
trino were obtained by considering effects of the emitted photons on the background 
radiation and big bang nucleosynthesis. Even today these effects lead to the most 
stringent constraint on long-lived radiatively decaying particles including massive 
neutrino. 

3.2.1. CMB spectral distortion 

As for the effect on the background radiation, the emitted photons (or electron- 
positron pairs) cause the spectral distortion of CMB as pointed out in Sato-Kobayashi 
1977. ^ If the neutrino decays early enough the photons injected in the decay are fully 
thermalized by photon number violating processes. In Ref.l), the bremsstrahlung 
(e+p — > e+p + j) was considered as the main photon number violating process. Af- 
terwards, double Compton scattering (e + 7 — > e + 27) was recognized to thermalize 
the background photons more efficiently by Lightman. 20 ) If we include the double 
Compton scattering, it is found that the injected photons are thermalized and the 
Planck distribution is recovered quickly for t < 10 6 sec. 

However, photons produced between decoupling of double Compton scattering 
(10 6 sec) and recombination (10 13 sec) may cause the spectral distortion in CMB. 
When the energy exchange between electrons and photons occurs sufficiently fre- 
quently by Compton scatterings (t < 10 8 sec), the spectrum of CMB photons be- 
comes Bose-Einstein distribution with chemical potential p. On the other hand, 
if the Compton scattering is less frequent (t > 10 9 sec) the spectrum distortion 
is described by y parameter.® The COBE satellite launched in 1989 showed that 
the CMB spectrum is that of a nearly perfect blackbody, from which the distortion 
parameters fi and y are stringently constrained as 21 ) 



H\ < 9 x 10~ 5 y < 1.2 X 10~ 6 . 



(3-9) 
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These distortion parameters are related to the energy injection Z\p 7 as 

Afr ^ { 0.714^ (10 6 sec < r„ h < 10 9 sec) 
Pi ~ \ % (10 9 sec < r Vh < 10 13 sec) ' 



(3-10) 



where Z\p 7 /p 7 is given by Eq. (|2-15p ( a factor (4/3) should be multiplied for decay 
into an electron-positron pair). Thus we can obtain the following constraint: 

( 10 3 sec(^) (m„<lkeV) 
r Vh <\ 10 6 sec (1 keV < m Uh < 30 MeV) , (3-11) 

1 10 12 sec(^) 4 (m Vh > 30 MeV) 

where we have taken account of "time dilation" , i.e. r e fj (effective lifetime)~ (Td/m Vh )T Uh 
in the small neutrino mass range. 22 )> 23 ) 

3.2.2. Big Bang Nucleosynthesis 

The BBN constraint discussed in Sato-Kobayashi 1977 was rather rough one 
and the more precise calculation was done by Miyama and one of the present author 
(KS), 22 ) which was later updated in Ref.24). In those studies, the abundances of 
the light elements (D, 3 He, 4 He, Li ) were estimated by calculation of the nuclear 
reaction network, and the result was compared with those inferred from observations. 
The resultant constraint is approximately given by 

Ty < / 3 >< 102 sec (^)" 2 S 2 MeV ) (342) 

Uh ~\ 10 2 sec (2 MeV < m Vh < 20 MeV) 

The above constraint is imposed mainly from entropy production due to the 
radiative decay of the massive neutrino. However, Lindley 25 ) pointed out that the 
radiative decay of the neutrino with mass larger than 10 MeV has another seri- 
ous effect on BBN. High energy photons with O(10) MeV produced in the decay 
can destroy the light elements which are synthesized in BBN through the following 
reactions: 

4 He + 7 -> T + p (19.8 MeV), (3-13) 
4 He + 7 -»• 3 He + n (20.6 MeV), (3-14) 
3 He + 7 -> B+p (5.5 MeV), (3-15) 
D + 7 -> n + p (2.2 MeV), (3-16) 



where the numbers in parentheses denote the threshold energies. High energy pho- 
tons and/or electrons injected into the cosmic plasma induce electro-magnetic show- 
ers and quickly thermalized, and how fast photons lose their energy crucially depends 
on whether the photon energy EL is high enough to scatter off the background pho- 
ton to produce electron-positron pairs (7 + 7bg e~ + e + ). For pair creation to 
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occur both energies of the background photon and the photon in the shower should 
be large, which leads to the condition, 

?r? 2 

E 7 >E th = ^. (3-17) 

Thus, high energy photons with EL > Eth can be thermalized by the photon-photon 
process so quickly that they have little chance to destroy the light elements. On the 
other hand, photons with lower energy cannot partake in the photon-photon process 
and are only thermalized by much slower Compton scatterings. Such photons can de- 
stroy the light element if their energy is above the thresholds of the photo-dissociation 
reactions. Using Eqs. (|3-13p - (13-17p it is found that when the heavy neutrino decays 
at t > 10 4 (10 6 ) sec, the photons in the induced electromagnetic showers can destroy 
D ( 4 He). This photo-dissociation effect was studied in details for radiative decay of 
massive neutrinos in Refs.24), 26), which leads to the following constraint: 

T^<10 4 sec (10 MeV < m„ ft < 1 GeV). (3-18) 

3.2.3. Constraint from Supernova 1987A 

So far we have reviewed the progress of the cosmological constraints on the ra- 
diatively decaying neutrino from the CMB distortion and BBN. Other astrophysical 
and experimental constraints are also discussed in the literature. Among these con- 
straints, here, we present the one from the supernova 1987A. SN1987A is a historic 
supernova since the neutrinos emitted its core were observed first in the human 
history. 

As well known, the neutrinos emitted from the supernova carry almost all ex- 
plosion energy ( ~ 10 53 erg) which is much larger than is visible (~ 19 erg). Thus, 
if even a small fraction of the neutrinos decay into photons (or electrons) inside the 
envelop of the progenitor ( < 3 x 10 12 cm ), the visible luminosity of the supernova is 
significantly changed, from which the stringent constraint on the lifetime is obtained 
as 27 ) 

f 10 7 sec(^) (m^ < 10 MeV) 

^ ~ 1 4 x 10 7 sec (^) 3/2 exp (-Jgfr) (m Uh > 10 MeV) 

One can see from Eqs. (J33D, (HES ), USED} , (l3~T8)) and (13091) that the 

radiative decay of the massive neutrino is almost forbidden. 

3.3. Non-radiative Decay 

In some models a massive neutrino can decay non-radiatively. For example, 
in the majoron (familon) model where a global lepton number (family) symme- 
try is spontaneously broken and associated Nambu-Goldstone boson called a ma- 
joron (familon) appears, a neutrino can decay into a lighter neutrino and a majoron 
(familon). Non-radiative decay was studied in Ref.29), and more recently in Ref.30). 
From these studies, it is known that non-radiative decay is also stringently con- 
strained by cosmology. 
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3.4. Number of Neutrinos 

The number of neutrino species N v was determined by measuring the decay 
width of the Z boson, which leads to 16 ) 

N v = 2.9840 ± 0.0082. (3-20) 

Therefore we now know that the only three types of neutrinos exist in nature. 

§4. Conclusion 

We have seen that the 1977 paper by Sato and Kobayashi 1 ) is considered as a 
historic paper which first showed that valuable information on neutrino properties 
can be obtained from cosmological consideration. The obtained constraints on ra- 
diative decay of a massive neutrino were updated by a number of authors and led 
to the important conclusion that the neutrino radiative decay is almost forbidden. 
The value of the paper is not limited in neutrino physics. The arguments also apply 
to more generic particles that decay on cosmological time scales as the subsequent 
studies showed. 

The paper initiated many studies to use cosmology as a laboratory of particle 
physics or to use particle physics to explore the very early universe. Today this 
research field is called particle cosmology and has achieved brilliant success. For 
example, the idea of inflation 31 )' 32 ) was born during effort to make monopole in 
the grand unified theories compatible with cosmology. Another example is interplay 
between supersymmetry and cosmology. Theories based on supersymmetry (SUSY) 
are considered as the most promising candidate for models beyond the standard 
particle physics, and predict many new particle some of which have long lifetimes. 
One of such long-lived particles is stable and account for the dark matter of the 
universe, 33 ^ while others often cause serious cosmological problems like gravitino 
problem. 34 ) In any way, it is now the standard procedure to test particle physics 
model by cosmology. 

In concluding this article, we would like to quote the final sentence of Ref.l), 

It is very interesting that valuable information on particle physics can be 
derived from cosmological arguments in spite of large uncertainties inher- 
ited there. 

K. Sato and M. Kobayashi (1977) 
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